Multi--dimensional Cosmological Radiative Transfer with a Variable Eddington Tensor Formalism

We present a new approach to numerically model continuum radiative transfer based on the Optically Thin Variable Eddington Tensor (OTVET) approximation. Our method insures the exact conservation of the photon number and flux (in the explicit formulation) and automatically switches from the optically thick to the optically thin regime. It scales as N logN with the number of hydrodynamic resolution elements and is independent of the number of sources of ionizing radiation (i.e. works equally fast for an arbitrary source function). We also describe an implementation of the algorithm in a Soften Lagrangian Hydrodynamic code (SLH) and a multi--frequency approach appropriate for hydrogen and helium continuum opacities. We present extensive tests of our method for single and multiple sources in homogeneous and inhomogeneous density distributions, as well as a realistic simulation of cosmological reionization.Comment: Accepted for publication in New Astronomy. Color GIF versions of figures 6, 7, 8, and 11 are available at http://casa.colorado.edu/~gnedin/PAPERPAGES/rt.htm

Distinctive rings in the 21 cm signal of the epoch of reionization

It is predicted that sources emitting UV radiation in the Lyman band during the epoch of reionization (EoR) showed a series of discontinuities in their Ly-alpha flux radial profile as a consequence of the thickness of the Lyman line series in the primeval intergalactic medium. Through unsaturated Wouthuysen-Field coupling, these spherical discontinuities are also present in the 21 cm emission of the neutral IGM. In this article, we study the effects these discontinuities have on the differential brightness temperature of the 21 cm signal of neutral hydrogen in a realistic setting including all other sources of fluctuations. We focus on the early phases of the EoR, and we address the question of the detectability by the planned Square Kilometre Array. Such a detection would be of great interest, because these structures could provide an unambiguous diagnostic for the cosmological origin of the signal remaining after the foreground cleaning procedure. Also, they could be used as a new type of standard rulers. We determine the differential brightness temperature of the 21 cm signal in the presence of inhomogeneous Wouthuysen-Field effect using simulations which include (hydro)dynamics and both ionizing and Lyman lines 3D radiative transfer with the code LICORICE. We find that the Lyman horizons are clearly visible on the maps and radial profiles around the first sources of our simulations, but for a limited time interval, typically \Delta z \approx 2 at z \sim 13. Stacking the profiles of the different sources of the simulation at a given redshift results in extending this interval to \Delta z \approx 4. When we take into account the implementation and design planned for the SKA (collecting area, sensitivity, resolution), we find that detection will be challenging. It may be possible with a 10 km diameter for the core, but will be difficult with the currently favored design of a 5 km core.Comment: 10 pages, 10 figures; v2: Section 5.5 rewritten; some new references added; accepted for publication in Astronomy and Astrophysic

Constrained semi-analytical models of Galactic outflows

We present semi-analytic models of galactic outflows, constrained by available observations on high redshift star formation and reionization. Galactic outflows are modeled in a manner akin to models of stellar wind blown bubbles. Large scale outflows can generically escape from low mass halos (M<10^9 M_sun) for a wide range of model parameters but not from high mass halos (M> 10^{11} M_sun). The gas phase metallicity of the outflow and within the galaxy are computed. Ionization states of different metal species are calculated and used to examine the detectability of metal lines from the outflows. The global influence of galactic outflows is also investigated. Models with only atomic cooled halos significantly fill the IGM at z~3 with metals (with -2.5>[Z/Z_sun]>-3.7), the actual extent depending on the efficiency of winds, the IMF, the fractional mass that goes through star formation and the reionization history of the universe. In these models, a large fraction of outflows at z~3 are supersonic, hot (T> 10^5 K) and have low density, making metal lines difficult to detect. They may also result in significant perturbations in the IGM gas on scales probed by the Lyman-alpha forest. On the contrary, models including molecular cooled halos with a normal mode of star formation can potentially volume fill the universe at z> 8 without drastic dynamic effects on the IGM, thereby setting up a possible metallicity floor (-4.0<[Z/Z_sun]<-3.6). Interestingly, molecular cooled halos with a ``top-heavy'' mode of star formation are not very successful in establishing the metallicity floor because of the additional radiative feedback, that they induce. (Abridged)Comment: 27 pages, 31 figures, 2 tables, pdflatex. Accepted for publication in MNRA

Simulations of Early Baryonic Structure Formation with Stream Velocity: I. Halo Abundance

It has been recently shown that the relative velocity between the dark matter and the baryons (vbc) at the time of recombination can affect the structure formation in the early universe (Tseliakhovich & Hirata 2010). We statistically quantify this effect using large cosmological simulations. We use three different high resolution sets of simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between dark matter and baryons. We show that the total number density of halos is suppressed by ~ 20% at z = 25 for vbc = 1{\sigma}(vbc), where {\sigma}(vbc) is the variance of the relative velocity, while for vbc = 3.4{\sigma}(vbc) the relative suppression at the same redshift reaches 50%, remaining at or above the 30% level all the way to z = 11. We also find high abundance of "empty halos", i.e., halos that have gas fraction below half of the cosmic mean baryonic fraction fb. Specifically we find that for vbc = 1{\sigma}(vbc) all halos below 10^5M\odot are empty at z \geq 19. The high abundance of empty halos results in significant delay in the formation of gas rich mini-halos and the first galaxies.Comment: 7 pages, 8 figures, accepted to Ap

The simulated 21 cm signal during the epoch of reionization : full modeling of the Ly-alpha pumping

The 21 cm emission of neutral hydrogen is the most promising probe of the epoch of reionization(EoR). In the next few years, the SKA pathfinders will provide statistical measurements of this signal. Numerical simulations predicting these observations are necessary to optimize the design of the instruments. The main difficulty is the computation of the spin temperature of neutral hydrogen which depends on the gas kinetic temperature and on the level of the local Lyman-alpha flux. A T_s >> T_cmb assumption is usual. However, this assumption does not apply early in the reionization history, or even later in the history as long as the sources of X-rays are too weak to heat the intergalactic medium significantly. This work presents the first EoR numerical simulations including, beside dynamics and ionizing continuum radiative transfer, a self-consistent treatment of the Ly-alpha radiative transfer. This allows us to compute the spin temperature more accurately. We use two different box sizes, 20 Mpc/h and 100 Mpc/h, and a star source model. Using the redshift dependence of average quantities, maps, and power spectra, we quantify the effect of using different assumptions to compute the spin temperature and the influence of the box size. The first effect comes from allowing for a signal in absorption. The magnitude of this effect depends on the amount of heating by hydrodynamic shocks and X-rays in the intergalactic medium(IGM). The second effects comes from using the real, local, Lyman-alpha flux. This effect is important for an average ionization fraction of less than 10%: it changes the overall amplitude of the 21 cm signal, and adds its own fluctuations to the power spectrum.Comment: 20 pages, 16 figures, 2 tables, To be publish A&A. High resolution version available at http://aramis.obspm.fr/~baek/21cm_Lya.pd

Discovery of a massive variable star with Z=Zo/36 in the galaxy DDO 68

The Local Volume dwarf galaxy DDO 68, from the spectroscopy of its two brightest HII regions (Knots 1 and 2) was designated as the second most metal-poor star-forming galaxy [12+log(O/H)=7.14]. In the repeated spectral observations in 2008 January with the 6-m telescope (BTA) of the HII region Knot 3 [having 12+log(O/H)=7.10+-0.06], we find a strong evidence of a transient event related to a massive star evolution. From the follow-up observation with the higher spectral resolution in 2008 February, we confirm this phenomenon, and give parameters of its emission-line spectrum comprising of Balmer HI and HeI lines. The luminosities of the strongest transient lines (Ha, Hb) are of a few 10^36 erg s^-1. We also detected an additional continuum component in the new spectrum of Knot 3, which displays the spectral energy distribution raising to ultraviolet. The estimate of the flux of this continuum leads us to its absolute V-band magnitude of ~-7.1. Based on the spectral properties of this transient component, we suggest that it is related to an evolved massive star of luminous blue variable type with Z=Zo/36. We briefly discuss observational constraints on parameters of this unique (in the aspect of the record low metallicity of the progenitor massive star) event and propose several lines of its study.Comment: 6 pages, 5 Postscript figures, to appear in MNRAS Letters in June 2008 issu

Modeling The Large Scale Bias of Neutral Hydrogen

We present new analytical estimates of the large scale bias of neutral hydrogen (HI). We use a simple, non-parametric model which monotonically relates the total mass of a halo M_tot with its HI mass M_HI at zero redshift; for earlier times we assume limiting models for the HI density evolution consistent with the data presently available, as well as two main scenarios for the evolution of our M_HI - M_tot relation. We find that both the linear and the first nonlinear bias terms exhibit a strong evolution with redshift, regardless of the specific limiting model assumed for the H I density over time. These analytical predictions are then shown to be consistent with measurements performed on the Millennium Simulation. Additionally, we show that this strong bias evolution does not sensibly affect the measurement of the HI power spectrum.Comment: 9 pages, 7 figures. Accepted by ApJ. New version: important changes from considering blue galaxy fraction, but conclusions remain the same. Fixed typo

The Mass Function of Dark Halos in Superclusters and Voids

A modification of the Press-Schechter theory allowing for presence of a background large-scale structure (LSS) - a supercluster or a void, is proposed. The LSS is accounted as the statistical constraints in form of linear functionals of the random overdensity field. The deviation of the background density within the LSS is interpreted in a pseudo-cosmological sense. Using the constraints formalism may help us to probe non-trivial spatial statistics of haloes, e.g. edge and shape effects on boundaries of the superclusters and voids. Parameters of the constraints are connected to features of the LSS: its mean overdensity, a spatial scale and a shape, and spatial momenta of higher orders. It is shown that presence of a non-virialized LSS can lead to an observable deviation of the mass function. This effect is exploited to build a procedure to recover parameters of the background perturbation from the observationally estimated mass function.Comment: 23 pages, 6 figures; to be appeared in Astronomy Reports, 2014, Vol. 58, No. 6, pp. 386-39

21-cm cosmology

Imaging the Universe during the first hundreds of millions of years remains one of the exciting challenges facing modern cosmology. Observations of the redshifted 21 cm line of atomic hydrogen offer the potential of opening a new window into this epoch. This would transform our understanding of the formation of the first stars and galaxies and of the thermal history of the Universe. A new generation of radio telescopes is being constructed for this purpose with the first results starting to trickle in. In this review, we detail the physics that governs the 21 cm signal and describe what might be learnt from upcoming observations. We also generalize our discussion to intensity mapping of other atomic and molecular lines.Comment: 64 pages, 20 figures, submitted to Reports on Progress in Physics, comments welcom

Visual Field Endpoints Based on Subgroups of Points May Be Useful in Glaucoma Clinical Trials: A Study With the Humphrey Field Analyzer and Compass Perimeter

PRECIS: Visual field endpoints based on average deviation of specific subsets of points rather than all points may offer a more homogenous dataset without necessarily worsening test-retest variability and so may be useful in clinical trials. PURPOSE: To characterize outcome measures encompassing particular subsets of visual field points and compare them as obtained with Humphrey (HVF) and Compass perimeters. METHODS: 30 patients with imaging-based glaucomatous neuropathy performed a pair of 24-2 tests with each of 2 perimeters. Non-weighted mean deviation (MD) was calculated for the whole field and separate vertical hemifields, and again after censoring of points with low sensitivity (MDc) and subsequently including only "abnormal" points with total deviation probability of <5% (MDc5%) or <2% (MDc2%). Test-retest variability was assessed using Bland-Altman 95% limits of agreement (95%LoA). RESULTS: For the whole field, using HVF, MD was -7.5±6.9▒dB, MDc -3.6±2.8▒dB, MDc5% -6.4±1.7▒dB and MDc2% -7.3±1.5▒dB. With Compass MD was -7.5±6.6, MDc -2.9±1.7▒dB, MDc5% -6.3±1.5, and MDC2% -7.9±1.6. The respective 95% LoA were 5.5, 5.3, 4.6 and 5.6 with HVF, and 4.8, 3.7, 7.1 and 7.1 with Compass. The respective number of eligible points were 52, 42±12, 20±11 and 15±9 with HVF, and 52, 41.2±12.6, 10±7 and 7±5 with Compass. With both machines, standard deviation (SD) and 95%LoA increased in hemifields compared to the total field, but this increase was mitigated after censoring. CONCLUSIONS: Restricting analysis to particular subsets of points of interest in the visual field after censoring points with low sensitivity, as compared with using the familiar total field mean deviation, can provide outcome measures with a broader range of mean deviation, a markedly reduced SD and therefore more homogenous dataset, without necessarily worsening test-retest variability